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The group Ea for grand unification is combined witfi the generation symmetry group SO(3)g. 
The coupling matrices in the Yukawa interaction are identified with the vacuum expectation values 
of scalar fields which are representations of the generation symmetry. These values determine the 
hierarchy of the fermions as well as their mixings and CP- violation. This generation mixing appears 
in conjunction with the mixing of the standard model fermions with the heavy fermions present 
in the lowest representation of i?6. A close connection between charged and neutral fermions is 
observed relating for instance the CKM mixings with the mass splittings of the light neutrinos. 
Numerical fits with only few parameters reproduce quantitatively all known fermion properties. 
The model predicts an inverted neutrino hierarchy and gives rather strict values for the light and 
heavy neutrino masses as well as for the 0^2(3 decay parameter. It also predicts that the masses of 
the two lightest of six 'right handed' neutrinos lie in the low TeV region. 
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I. INTRODUCTION 

The origin of the properties of quarks and leptons, the 
masses and the mixings of the three generations, forms 
still an open problem in particle physics. Grand unified 
theories [l| provide a general understanding of the struc- 
ture and the quantum numbers of the standard model 
states and suggest new ideas. The hope is to find a 
consistent scheme which provides for intimate relations 
between the many observables. Lagrangians considered 
for a single generation have the fundamental property 
of chiral symmetry, i.e. they are invariant under chi- 
ral transformations of the fermion fields taken together 
with appropriate transformations of the Higgs fields. For 
more generations chiral symmetry requires an extended 
symmetry which relates the different generations: the 
fermions should be members of a non-Abelian genera- 
tion symmetry.. There are two possibilities to enforce 
generation symmetry: one can either enlarge the number 
of Higgs fields by giving them generation indices or one 
can identify the coupling matrices in front of the Higgs 
fields as vacuum expectation values of new scalar fields 
carrying generation quantum numbers. 

We choose here the second alternative in connection 
with the grand unified symmetry group Eq 0-01 ■ For 
three generations the coupling matrices are 3x3 matri- 
ces. The corresponding 9 scalar fields are taken to be 
hermitian and are expressed in terms of the hermitian 
3x3 matrix <I>(x) 

^x) = x{x)+iax) . (LI) 



Here x denotes a symmetric and ^ an antisymmetric 3x3 
matrix. Clearly, this choice implies that conventional 
Yukawa interactions containing these fields are effective 
ones with dimension 5 and thus have to be understood 
on a deeper level. This can be done by the introduction 
of additional heavy spinor fields. We will see, that by 
integrating out these heavy fields one finds interesting 
consequences for the relation between quark and neutrino 
mass matrices. 

The introduction of the hermitian matrix field <f>(x) 
coupled to the fermion fields suggests to use the group 
50(3)g to describe the generation symmetry. In addition 
we will make use of a discrete parity like symmetry, gener- 
ation parity Pg. From the point of view of chiral symme- 
try the use of SU{i)g instead of SU{2)g ~ SO{3)g would 
be more consequent. We will not treat this extension 
because the subgroup S0{3)g of SU{3)g, in combination 
with Pg, is sufficient to reach our aims. Together with 
the unification group Eq it can be dealt with in a very 
economical way. In the literature non-Abelian continu- 
ous 01 and discrete Q fiavor symmetries have been dis- 
cussed and applied in various models. In our approach all 
fermion fields are taken to transform as 3-vectors under 
the generation group SO{3)g. Consequently, the sym- 
metric part of $ has to transform as 1 -I- 5 and the an- 
tisymmetric part as 3. Spontaneous symmetry breaking 
of this generation symmetry leads to vacuum expectation 
values of $. By an orthogonal transformation the sym- 
metric matrix (x) can be taken diagonal. As we will see 
its elements describe the up quark hierarchy: 
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Here M denotes the scale at which the appropriate ef- 
fective Yukawa interaction of dimension 5 is formed. We 
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take the mass ratios rric/rnt = m„/TOc = to be vahd 
at a high scale. At the weak scale Mz these ratios are 
somewhat modified. Taking a ~ 0.050 gives good agree- 
ment with the experimental mass determinations. In the 
following we will use this parameter for expressing small 
quantities. 

There remains only the vacuum expectation value of 
the antisymmetric matrix to produce all mixing and 
CP violating properties of quarks and leptons. (The use 
of a purely antisymmetric and hermitian mixing matrix 
was suggested in Q). We take for ^ the particular form 
which was abstracted from the analysis of the fermion 
spectrum and the CKM matrix in Q using Eq symmetry 
for grand unification. It has a particular symmetry: it 
changes sign by exchanging the second with the third 
generation: 

1 / a -a \ 

— (iO =A = i \ o 1/2 . (1.3) 

\ cr -1/2 / 

The relative factor between the (1,2) and (2,3) elements 
determines the particle mixings. Its value 2a gives a 
good description of the mixings of quarks and neutrinos 
as will be shown in this article. The coupling constant 
\a can be incorporated into the mass scale M' (except 
for renormalization group considerations). 

As already mentioned above, besides the SO{?))g sym- 
metry we also introduce a parity symmetry Pg, 'gener- 
ation parity'. The (extended) standard model fermions 
are taken to be even under this symmetry while the Higgs 
fields directly coupled to these fermion fields as well as 
the fields x and ^ have negative generation parity. To ob- 
tain the spontaneous symmetry breaking of 5*0(3)^ x Pg 
we use in the Lagrangian for the field $ SO{i)g x Pg in- 
variant potentials up to 4th order in the fields x ^-i^d 
^. Adding also the SO{i)g x Pg invariant Coleman- 
Weinberg potential a complete breaking of the genera- 
tion symmetry can be achieved. Moreover, by selecting 
properly the coefficients of these potentials the numerical 
values of (|L2p and (jl.3p give the absolute minimum of the 
total potential. 

It is clearly a challenge to obtain the masses, mixings 
and CP properties of all fermions with only the two gen- 
eration matrices G and A together with a few vacuum ex- 
pectation values of the corresponding Higgs fields. Suited 
for this task is the grand unification symmetry Eq Q. 
Here we can postulate that particle mixings are caused 
by the mixing of the standard model particles with the 
heavy particles occurring in the lowest representation 27 
of Eq. The up quarks cannot mix since they have no 
heavy partners in this representation. Thus, it follows 
immediately that - apart from a constant factor - their 
mass spectrum is simply given by (x) /M = G [as pre- 
sented in ((L2)) ]. 

As we will see the breaking of generation symmetry 
and Eq symmetry gives a hold on the complete fermion 
spectrum which includes new heavy fermions. In partic- 
ular, we find a very strong hierarchy for the heavy right 



handed neutrinos: their mass matrix is proportional to 
the square of the up quark mass matrix. By integrating 
out these heavy neutrinos, their masses appear in the 
denominator compensating the square of the Dirac mass 
matrix in the nominator. This mechanism gives the light 
neutrino spectrum a less pronounced hierarchy than the 
one of the up quarks. A possible consequence is an in- 
verted neutrino spectrum. Such an interesting situation 
is not obtained in more conventional treatments of grand 
unified theories. It will be seen that very few parame- 
ters are sufficient to describe all the known properties of 
charged and neutral fermions in a quantitative way. 

The use of Eq as a grand unified theory has many 
virtues The fermions arc in the lowest repre- 

sentation of this group and an elegant cyclic symmetry 
connects quark fields, Icpton fields and anti-quark fields. 
In particular, as shown in Q, Eq combines the mix- 
ings among fermion generations with the mixing of the 
standard model particles with heavy charged and neutral 
states. 

Thus, our starting symmetry at the GUT scale is 

Eq X S0i3)g X Pg . (1.4) 

In order to introduce notations and conventions we add 
at this place a short description of the 27 representation 
of Eq. In the Eq grand unification model the fermions 
are contained in the 27 representation of the group, i e. 
they are described by 27 two component (left handed) 
Weyl fields for each generation: 

C , « = 1,2,3 , r = l,... ,27. (1.5) 

r denotes the Eq fiavor index and a labels the genera- 
tions. These fields - even under Pg - describe the fermions 
of the standard model plus additional quark and anti- 
quark fields with the same charge as the down quarks and 
new heavy charged and neutral leptons. All fermions are 
in singlet and triplet SU{3) representations of the maxi- 
mal subgroup of Eq 

SU{3)l X SU{3)r X SUi3)c = G333 , (1.6) 

which plays an important role in our approach. In terms 
of G333 we have 

27 = Ql{x) + L{x) + Qr{x) , (1.7) 

where the quantum number assignments are: 

Quarks : Ql{x) = (3, 1, 3) , 
Leptons : L{x) = (3,3,1) , 
Anti- quarks: Qr{x) ^ (1,3,3) . (1.8) 

For each generation one has 

fu^X / L\ E- e- \ 

[QlT^ = M° , Ll = \ E+ Ll V \ , 

Ve+ Llj 

{Qnt = («a, 4, Da) , (1.9) 
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where i,k,a = 1,2,3. In this description 5C/(3)l acts 
vertically (index i) and SU{3)ii horizontally (index k) 
and a is a color index. It is seen, that in each generation 
there exist 12 new fields extending the standard model: 
colored quarks D,D, charged leptons E^, and 4 new 
neutral leptons v,L^,L\,L^, which all must correspond 
to heavy, but not necessarily very heavy, particles. 

The group G333 can serve as an intermediate gauge 
symmetry below the Eq breaking scale (the unification 
scale). It can be unbroken only at and above the scale 
where the two electro- weak gauge couplings combine. In 
the non-supersymmetric Eq model, which we adopt here, 
this point occurs at M/ ~ 1.3 • 10^'^ GeV according to 
the extrapolation of the standard model couplings. In- 
terestingly, this is just the scale relevant for the neu- 
trino masses using the see-saw mechanism. Above this 
scale the united electroweak couplings and the QCD cou- 
pling run at first separately until the electroweak cou- 
pling bends to meet the QCD coupling at the Eq unifi- 
cation point. This happens at the scale « 3 • 10^^ GeV 
(see Fig. [1] and more details in Appendix lA ip . 

Before Eq symmetry breaking equivalent forms of (|I.9p 
can be obtained by applying left and right SU{i) lA-spm 
rotations. We fix the basis by using vacuum expectation 
values (VEV) of the lowest Higgs representation, the 27 
of Eq. In this representation possible vacuum expecta- 
tion values are restricted to the 5 neutral members sitting 
in positions corresponding to the ones of the neutral lep- 
tons in p.9|) . Below the scale Mi only two light SU{2)l 
doublets arc assumed to be active. We incorporate them 
in two different scalar 27-plet fields. Thus we introduce 
two scalar fields and with the following trans- 
formation properties under E^ x 50(3)^ x Pg 

(27,1,-) , 7?27^ (27,1,+) . (1.10) 

Only can couple to the ifj fields in the form x ij'^Hip. 
Thus, it is convenient to choose a basis in which the VEV 
{H27)\. forms a diagonal matrix. In the scalar potential 
of the Eq Lagrangian the bilinear terms with respect to 
both scalar 27- plets can contain the singlet part X(i) of 
the fields %: 

fi^Hl,H27 + 11^HIH27 + Ai'X(l) (^17^27 + ^^27^17) ■ 

(1.11) 

After develops a VEV, and for fi' sufficiently small 
(i.e. Im^~A*^I)i the mass eigcnstates i/" and 

if following from p.lip are linearly related to H27 and 
H27 

H27 = i/" + zW^ , and H27 = -zff " + H"^ , 

with z = -J^^p^ , \z\ « 1. (1.12) 

(m - M ) 

We take the hght up type doublet to be in i/": 
((ff")J,(ff")^) and the light down type doublet 
{{H'^)l, {H'^)l) to occur in H'^. These Higgs doublets are 
important for generating the correct fermion mass pat- 
tern. ((i?27)i) — ~ determines the scale of 



the up quarks and ((^27)2) — ^{{H'^)'^) = ze^ the scale 
of the down quarks and charged leptons of the standard 
model. £2 is expected to be of the same order of mag- 
nitude as e\. The factor z, which vanishes before Pg- 
symmetry breaking, is responsible for the small values of 
bottom and tau masses compared to the top mass. Com- 
ponents of 77"", H'^ which arc standard model singlets 
have large VEVs. For instance ((7^27)3) - ((^f")!) = 4 
provide high (Dirac type) masses for all new quarks and 
leptons with the exception of the 'right handed' neutri- 
nos D = L2 and L3. Their masses arise from a different 
mechanism involving both H27 and 7?27- 63 is not a new 
scale parameter but is determined by the onset of G333, 
the meeting point of the electroweak gauge couplings: 
gi{Mj) = g2{Mi) with the result (see Appendix I A ip : 

Mi = 1.27-10" GeV , = , = 2.27-10" GeV . 

.92 (M7) 

(1.13) 

We will see that this value of eg provides for the correct 
mass scale of the light neutrinos. 

Due to the linear Yukawa coupling all fermions - ex- 
cept the 'right handed' neutrinos L2 and - arc Dirac 
particles and (so far) have the same G hierarchy as the 
up quarks. The hierarchy of the 'right handed' neutri- 
nos, on the other hand, becomes super strong due to the 
combined action of H27 and 7?27- Particle mixings will 
occur by taking into account the antisymmetric genera- 
tion matrix A combined with an antisymmetric (in Eq 
indices) Higgs field. These mixings modify the hierarchy 
of down quarks and charged and neutral leptons in agree- 
ment with the experimentally observed particle spectra. 

II. GENERATION SYMMETRY 

For the construction of the Yukawa interaction sym- 
metric under Eg x SO{3)g x Pg the existence of a set of 
new fields is necessary. First of all, as mentioned in the 
introduction, we introduce scalar fields represented by 
the 3x3 matrix $(2:). The members of the symmetric 
part X transforms as (1, 1-1-5, — ) [under Eq x SO{3)g x Pg] 
and the ones of the antisymmetric part ^ according to 
(1,3,—). To have a renormalizable interaction we need 
new heavy vector like fermionic fields F,F and F',F' 
which transform as 

F^(27,3,-), F^(27,3,-), 
F'^(27,3,-), F'~(27,3,-). (11.14) 

The Yukawa interaction involving the Higgs fields 
arises from the vertices 

(V'^(5(i)X(i)+5(5)X(5))i^) , M{F^F), {F^H27^) ■ 

(11.15) 

In pi.lSp we used matrix notation with regard to gen- 
eration indices, but suppressed Eq indices and Clebsch 
Gordan coefficients. and ^(5) denote the singlet and 
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5-plct parts of Xj respectively. By integrating out the 
massive fields F, F one gets the wanted effective interac- 
tion 



1 

M 



(V'^(5(l)X(l) + 9{5)X(5))H27^) 



(11.16) 



The corresponding diagram is shown in Fig. [2^. Clearly, 
the antisymmetric matrix ^ does not contribute to the 
symmetric (in Eg indices) field H27 ■ It couples, however 
to the antisymmetric Higgs representation iJ^(351) of Eg 
with negative generation parity. The vertices are 



{{i^^9i3)^(F') , M'iF'^F') , (F'^i/^35i)^) 



(11.17) 



The corresponding diagram containing the heavy fields 
F',F' is shown in Fig. Eb. 

The Pg symmetry forbids couplings such as ipHip , 
FHip, FFH etc. In the following we will not separate 
the singlet and 5-plet parts of x but simply use the com- 
bination X 8-8 occurring in $ (i.c g^i^ = 1/3 5(5) ) and 
absorb (7(5) and the coupling 3(3) in the values of M and 
Af, respectively. Now we can introduce number valued 
generation matrices by taking the vacuum expectation 
values of x C described in the introduction: 



M 



{X)=G 



1 



(11.18) 



As mentioned before, the generation matrix G in Eq. 
pL18| combined with vacuum expectation values of the 
Higgs field gives Dirac masses to all fermions ex- 
cept the two heavy leptons L2 and L3. The latter re- 
quire Higgs fields transforming as (6, 6) with respect to 
SU{3)l X SU{3)fi, components, which are not contained 
in H and Ha- Instead of introducing another high di- 
mensional Higgs field it is plausible to use the Higgs field 
H27 together with H27 ■ This avoids the appearance of 
a new unknown generation matrix. With the help of the 
generation symmetry it is possible to derive the relevant 
generation matrix of the 'composite' (6, 6) plet in terms 
of the G matrix. This increases the predictive power of 
the model. To obtain the corresponding effective interac- 
tion from a renormalizable interaction, another massive 
Dirac field is required. It is a vector in generation space 
and an Eg singlet: iV(l,3,+), iV(l,3,-) with N even 
under Pg and N odd under Pg. In addition we need a 
total singlet field (1,1,—) which is odd under Pg. It can 
be identified with the singlet part X(i) of x with vacuum 
expectation value = A/jv ~ M. The vertices are 



f^hIjN 



( nhIjF 



{X(i)){NN)^Mn{NN), 

(11.19) 

Integrating out the fields N, N the vertex 
(F'^ H^){H^ F)/Mn emerges. Below the masses of 
the F states the effective Yukawa interaction 



1 1 

M2 



r xHl7] {H^7Xi^ 



(H.20) 



is generated. The corresponding diagram is shown in 
Fig. [3l Because in (|II.20p the Eq indices are not shown, 
one has to keep in mind that the combinations i/jH^ and 
tpW are Eq singlets. It is easily seen that by this inter- 
action only neutral leptons can get masses, notably the 
right handed neutrinos L2 and L3 . They are coupled to 
the large elements of H27 and H27, which are standard 
model singlet fields and thus will get large vacuum ex- 
pectation values. The most interesting feature is however 
the appearance of the square of x/M . It implies that af- 
ter generation symmetry breaking the mass hierarchy of 
these heavy neutrinos is dramatic, namely equal to G^, 
the square of the mass hierarchy of the up quarks which 
is already a very strong one. 

i?27 and J which both appear now in the Yukawa 
interaction, can be expressed in term of the mass eigen- 
statcs H"- and H"^. From (|LT2|1 we have H27 = H^ + zH"^ 
and H27 = -zH" + H<^, with the mixing parameter 
|z| < 1 . We choose the VEVs of if" and H'^ as fol- 
lows 
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(11.21) 



Thus, the VEV of the Higgs field H27 has the diagonal 
form taken by convention 



H27) = (i?" 








(H.22) 



while i?27's VEV has the structure 



(Ho 



(H^-zHn 



-ze 









(11.23) 



Since (H27) has no off diagonal element one expects a 
large (3,2) element for H27. We take ~ Mn ^ M, 
where M is identified with 63 ~ A//. In principle, how- 
ever, £2 together with M and A/jy could be of a lower 
scale (but still 3> Mz). The dominant VEVs of H27 and 
H27 fix now the mass terms of the 'right handed' neutri- 
nos L2 and L3 

F{^^^}G^((LlfLl) + F^^^^^G^((LlfLl) , 



with F{2.3} ~ ^ 



(H.24) 

^{3,3} jg ^Yic only Majorana mass term occurring so far. 

The effective Yukawa interaction below (M, AIn , M') 
now reads 



^{3,3} 



-z(eir 



-ze 



G^li {i,'-^H27^^) + (^"^i/A35i^'') + 

k-N (^').. [r^iil^iiU^) . (H.25) 
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This effective Yukawa interaction together with the VEV 
configurations (jll.22p and pL23p contains all the neces- 
sary information about the generation structure: Gen- 
eration hierarchy and generation mixing are now com- 
pletely fixed. At this stage we do not need to specify the 
scales of M, Mm — M, M' ^ Mz- However, these scales 
become important when we take renormalization effects 
into account. 

Let us now discuss the breaking of the generation sym- 
metry 50(3)3. The Lagrangian for the field $ is 



U = 2Tr {(a^$ _ c[B^, $]) (a^a> - c[b^, m-yix,0 , 



with Bf^ = BJ^t^, where Bj^ denote the vector potential 
and the 3 antisymmetric generators of the generation 
symmetry SO{3)g. V{xtO stands for an S0{3)g invari- 
ant potential. Gauge invariance allows to choose X: the 
symmetric part of to be a diagonal matrix with 3 real 
elements. This defines a direction in symmetry space for 
a possible spontaneous symmetry breaking. Normalizing 
$ for this basis one can write 



tential 



XI 
X2 
X3 




(11.26) 



In this form the vertices we have shown above simplify 
drastically since the orthogonal transformation diagonal- 
izing X can be absorbed by the fields i/' , F and F' . 



The scalar potential V{x,£,) in pi.26[) can easily be 
made to have a minimum for specific values of the three 
invariants: the trace of x^, the trace of the square of the 
traceless part of x and the trace of At this mini- 
mum we have then two relations for the 3 fields form- 
ing X and one relation for the other 3 fields forming ^. 
There still remains the freedom of such SO{3)g transfor- 
mations which do respect this minimum and keep x diag- 
onal. This remaining symmetry is necessarily a discrete 
subgroup of S0{3)g. It is the discrete symmetry group 
A4. As a generation symmetry has been suggested in 
many publications Q. In our approach this symmetry 
occurs naturally in connection with the starting symme- 
try S0{3)g xPg. It is the remaining symmetry after using 
appropriate potentials invariant under SO{3)g x Pg and 
the choice of a symmetry breaking direction. 

Radiative corrections add to the potential V{x, new 
S0{3)g invariant parts which are of logarithmic type. 
The Coleman- Weinberg potential is of this form. By 
including it the total potential can lead to a complete 
spontaneous symmetry breaking of SO{3)g which results 
in vacuum expectation values for x and ^ of the form 
(jl.2l) . p.3p . We demonstrate this here by using the po- 



v{x,0 



Tr[x' 



ci(Tr[x])^ 
2n2 



2 "''^ ' 4 
CT(Tr[x|])2 +C3Tr[x3] Tr[x] + (Tr[(*0 

CxCoTr[x.(iO^] Tr[x] + CxeiTr[x.i?.X-iC] 



647r2 



M|i In di 



M 



,t2 

Bl 

M2 



M|,ln(d2^ 



+Ml,ln[d3|f 



(11.27) 



Here xt denotes the traceless part of x and Msi stands 
for the vector boson masses as field dependent functions 
(we take in pl.27[l e = 1 for simplicity). The coefficients 
in pL27P can be tuned such that the SO{3)g x Pg sym- 
metry breaks spontaneously and produces G and A with 
a ~ 0.050. To achieve this one has to use the six rela- 
tions following from the first derivatives of the potential 
at the proposed minimum: (x) = MG and (^) = M'A. 
Because of the large hierarchy a high accuracy of this 
calculations is necessary. We use here M'/M = 10^ i.e., 
the mass scale for forming the antisymmetric matrix A 
is large compared to the one for forming the symmetric 
matrix G (this is required for the gauge coupling unifi- 
cation, see the renormalization group treatment in the 
appendix). Setting then e.g. di = 1, ^2 = 237854/10^ 
and da = 237850/10^ all other coefficients are determined 
by putting the first derivatives of the potential to zero. 
These coefficients are sufficiently small (ci,ct,C3,c^ < 1 
for instance) to allow perturbative treatments. The log 
terms are small near the minimum. With this choice of 
the potential the six eigenvalues of the 6x6 matrix for 
the second derivatives of the potential are all positive 
at the wanted values of ($). The minimum obtained is 
an absolute one (but degenerate with respect to different 
signs of the three elements in (^)). 

The above formulae for the vacuum expectation values 
and the corresponding potential would have a different 
form if we would have used other values for the cou- 
plings 5(1) and 3(5). However, also for this more general 
situation potentials can be constructed to produce the 
required spontaneous symmetry breaking. 

For the three vector boson masses at the minimum of 
the potential one gets 



1.40 eM 
0.505 eM' 
0.505 eM' 



(H.28) 



III. 



CHARGED FERMION MASSES AND 
MIXINGS 



As we have already mentioned, the diagonal generation 
matrix (x) (i.e. G) is taken such that the up quark masses 
have their observed hierarchy. The top quark mass at 
Mz is determined by the vacuum expectation value e\ 
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and the coupling At at this scale 

mt = e\\t = e\ = 170.9 GeV 



(III.29) 



The notation (the superscript uii of G) indicates that 
at low energy one has to distinguish between different 
channels. The renormalization group effects are treated 
in the Appendix. We take for Xt{Mz) a value such that 
the top coupling constant near Mi is close to one. 



XtiMz) = 1.30 



el 



131.5 GcV 



(III.30) 



Having for all diagonal G-couplings the same spacing at 
the scale Mj , we are able to calculate the spacing for each 
channel at the scale Mz- As obtained in the Appendix 
we find 

G""(Afz) = Diag (rr , 7- , l) • Xt{Mz) , 

with rt = 0.717 . (III.31) 

Comparing with the measured values of up and charmed 
quark masses suggest 



a = 0.050 . 



(III.32) 



Eq. (jIII.3ip gives for this value of cr iriu = 1.5 MeV and 
rric = 0.60 GeV in good agreement with experiment. We 
will use pil.32|l in all further calculations. 

The case for the down quark and charged lepton masses 
is not as straightforward. Here the light fermions will mix 
with the heavy states D and L. We take this as the source 
of the particle mixings. It involves the indices 2 and 
3 {U-spin mixing) of SU{3)l and SU{3)ji of the Higgs 
field Ha- According to Eq. pi.25p this mixing of light 
and heavy fermions occurs together with the generation 
mixings by the coupling matrix z(^) (i.e. the matrix A) . 

i) The dovifn quark masses and mixings 

The matrix elements for down quarks are part of a 6 x 6 
mass matrix. It has the form 



d,D 



d 

D 



P2A, 



D 



(IH.33) 



In this equation the constants /j stand for the vacuum 
expectation values of 77^(3, 3, 1): 



/; = ((H^)}), z,j = 2,3 



(IH.34) 



e| and /| have values of the order of the weak scale, 
while 63 and /I will be very large because they arise from 
standard model singlets. Further below we will find that 
the vacuum expectation values of the standard model 
singlet components of Hj\ , like /| , have values which are 
small compared to 63 (smaller than <J^e\). This allows 
to make use of the see-saw formula and also justifies the 



neglection of /| in pil.33p . Integration of the heavy D- 
states leads to the 3x3 down quark mass matrix at the 
scale Mz 



m4Mz) = zelG''^ + A''''fi 



(IH.35) 

We distinguished different G and A matrices when cou- 
pled to different channels {dd, dD) etc. Different renor- 
malization factors arise when starting from the original 
G and A matrices at high scales. The calculations are 
deferred to the appendices. The matrices appearing in 
PII.35P are defined in (f02|l and (f04)l . 

For the dd matrix element in pil.35[) with the genera- 
tion index (3,3) we write 



(IH.36) 



The down quark mass matrix leads to 7 obscrvablcs : 3 
mass eigenvalues, 3 mixing angles and the CP violating 
phase. According to pil.35[) we have 3 parameters for 
a fit of the experimental results. We use the notation 
of the various coupling constants provided in (|A.12p and 
(f04)) . Taking then 



ml = 2.95 GcV 



-0.23 GeV 



i-3 J-2 \dD\Dd 

^^i^ ^ ^ = 1.62 • 10-* GeV , (IH.37) 

one gets 

TO<i(Mz) =i 2.6 MeV , m^(M2) ~ 50 MeV , 
rubiMz) 2.89 GeV , 
IK.I ^ 0.228 , \Vcb\ =^ 0.042 , ^ 0.0039 . (IH.38) 
The angles in the quark unitarity triangle have the values 
aq ~ 97° , Pq ~ 23° , 7q ~ 60° . (IH.39) 

Our results are very satisfying. Quark masses, the CKM 
mixing angles and the three angles of the unitarity trian- 
gle have values within experimental errors [Til . [T3 | . 

ii) The charged lepton masses and mixings 

The charged lepton sector is constructed in a similar 
way. The light leptons mix with the heavy L's through 
the vacuum expectation values of the 7?yi(3, 6, 1) mul- 
tiplet (the vacuum expectation values of Ha{G,3, 1) are 
considered to be negligible because of its 6 representation 
for the SU{3)l symmetry): 



(H 



f 



(IH.40) 



With the abbreviations 

^2(1.3} ^^2^ ^2(1,2} ^ ^^2^ fn,3}^gl (111.41) 

and applying again the effective Lagrangian pi.25P the 
6x6 mass matrix for charged leptons has the form 

e+ E+ 
e- (-zelG-glA, g^A \ (IH.42) 



M, 



e.E 



E- 



4G 
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Integrating out the heavy L states leads at the scale fi ' 
Mz to the 3x3 matrix 



^ (III.43) 
Here appear again renormalization group coefficients de- 
fined in f02)) . ((XT4)) . We choose 



zGf^^ el ~ elK = = 1.593 GeV 
X^A^^gl ^ -0.1898 GeV , 
A 



9392 ■^A ■^A 



3.16- 10"^ GeV . (III.44) 



These three numbers determine the charged lepton 
masses and, as in the quark case, determine also the 3 
mixing angles and the CP violating phase. The charged 
lepton mixings are not directly observable, but will play 
an important role in the discussion of the neutrino prop- 
erties. The masses come out precisely 

me{Mz) = 0.487 MeV , m^{Mz) = 102.8 MeV , 

mr{Mz) = 1.747 GeV . (11145) 

For the mixing angles and the CP violating phase 7*^ we 
obtain 



'12 



^23 



5.1° 



^13 



0.35° 



7 



- -23° . 
(in.46) 

We have defined these four parameters in complete anal- 
ogy to the quark CKM mixing angles and the unitarity 
angle gamma. 

iii) Estimates of /|, /|, and g| 

So far only the combinations Ai-e2i /I /I ^'^^ 53 52 ^^'^ 
known numerically. However, there is one relation due to 
the known mass of the vector boson W of the standard 
model. It connects all vacuum expectation values which 
belong to SU{2)l doublets: 

{e\r+{4f+{fi)'+i9'2r+{fir+{9ir^m.iGcvr . 

(IH.47) 

Because of the nearly identical interaction of down 
quarks and charged leptons we can estimate the ratio 
from the fit values given in piL37p and (|III.44|) by 
taking the large standard model singlet VEV's /| and g| 

to be equal. One gets ^ ~ 0.5 which is close to m^/ml 

as one could have expected. Using Eq. piL47p the only 
parameter needed then is the ratio 3|/e| . This is a ratio 
of two weak scale quantities with the same weak isospin 
quantum numbers. We therefore take it to be equal to 
one. This choice gives a value for gf/^i small enough to 
justify the application of the see-saw mechanism we used 
above. It also provides for small values for the couplings 
At- and Af, necessary for renormalization group stability 
of the neutrino sector to be discussed in the appendices. 



All couplings and VEV's introduced are now fixed: 

fl ~ 39 GcV , gl = e2 76 GeV , 
zGf/ c^XrC^ hOGcV ^ 2 . 10-2 ^ ^ 4 . 10-2 



/I gl 



4.16 • 10" 



(HI. 



--3 "^3 



It is useful to set 



givmg 



/I .9i 



^3 '^3 



Xg ~ 0.033 . 



(IH.49) 



The factor cr^ removes the singularities in the mass ma- 
trices with respect to the formal limit cr — > . 



IV. NEUTRINO MASSES AND MIXINGS 

In each generation one has to deal with 5 neutral lep- 
tons (see p.9[) ). Thus, the matrix for neutral leptons is 
a 15 X 15 matrix. 



Ll 
L\ 



^2 



Ll 




Ll 


L\ 





~e\G 





-glA 


-e}G 





Ml 








Ml 


M2 





-^9lA^ 

















e\G 


M^ 



^2 




e\G 
Mo 





(IV.50) 

Each entry stands for a 3 x 3 matrix. In the 12 x 12 sub- 
matrix for the heavy leptons we neglected small terms 
like and g|. They play no role in the evaluation of the 
light neutrino properties. The following abbreviations 
are used: 



Mo ^elG 



Mi^ f'^'^'^^G^+FaA, M2=F^3'3>G^ 

(IV.51) 

Ml and M2 have the superstrong hierarchy G^ according 
to Cf and Eq. pL24)l . The only new element we had 
to introduce arises from the vacuum expectation value 
of _ff^(6, 3, 1) which wc argued to be negligible in the 
previous section. But here it appears together with very 
small elements occurring in G^ . Thus, the constant Fa 
defined according to 

Fa = {{Ha) {?,,?,}, 1) = /{3,3}i 

is expected to be tiny compared to F^^.s} gyen though it 
is a singlet with respect to standard model transforma- 
tions. Nevertheless it has to be kept as a parameter. 
It is useful to rewrite the matrix pV.50[l in the form 



Mr 



o 

M 



(IV.52) 
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where stands for the 3x3 zero block matrix, while 
and M are 3x12 and 12 x 12 matrices respectively. They 
read 



n = {-elG, 0, -.92 A 0) 



(IV.53) 



For consistency the quantity Xg, at this place responsible 
for the correction to the see-saw formula, must be small 
compared to 1 which is indeed the case as we have seen 
before. Finally, taking all together, one obtains for the 
light neutrino mass matrix 



M 



/ Ml \ 

A'h e\G 

A/o 

y e\G / 



(IV. 54) 



The matrix M contains the masses of the heavy neutral 
states which should be integrated out. This can be done 
analytically since gf/sl is sufficiently small according to 
pil.48p . In doing so the light neutrino 3x3 mass matrix 
is given by 

= -U'^mr^Q^U . (IV.55) 

In this expression the matrix U provides for a first order 
correction to the generalized see-saw result Q,M~^VL^ . 
One finds with sufficient accuracy 

1 



U ~ 1 



2nT 



(IV.56) 



The inverse of the matrix M is 



/ 1 Af 1 1 



M- 



1 
Ml 



M, 



1 ri i 

■T <-r M 



Ml 



1 



v 











1 

Mo 



From pV.53jl and (|IV.57j) one obtains 

mr^u^ = -(e 



/ 
(IV.57) 



A? 



Mf Ml 



Mf A'/o 



Ml 



(IV.58) 



Fig. 0] shows how the terms in pV.58p are generated. 
The generation of the first term in (jIV.58p is exhibited by 
the diagram of Fig. [4^, while the diagram of Fig. \^ and 
its transpose shows the formation of the second and third 
terms of pV.58p . respectively. In linear approximation 
with respect to Fa the matrix M^"^ becomes 



A/; 



F{2.3} 



1 

G2 



^{2,3} Q2 Q2 



(IV.59) 



The expression of lA appearing in Eq. (jIV.55|l can be 
approximated by 



U ~ 1 



-3,1 1 



1-^ 



(IV.60) 



^{2,3}" l^i;^{2,3} 

F{2,3} er G2 G ^ G G2 ■ 



1 



el ^^G + G 



A' 



U . (IV.61) 



In order to compensate the powers of a in the denomi- 
nators of (jIV.61|l we extract powers of a from Fa as we 
did for 52 ■ 



9l 



a Xo 



Fa 

i^{2,3} 



a^XA 



(IV.62) 



Like Xg also xa needs to be small compared to one. 

For a convenient description and discussion we use now 
the abbreviations 



X 



5 F{2,3} 



i;-{3,3} 

XgFi^'^} 



(IV.63) 



Let us at first discuss the case with only linear terms 
in Xg and cr, putting also xa ~ . The neutrino mass 
matrix takes then the very simple form 



p m 1 



— 1 1 
-i? 
-if 0^ 



(IV. 64) 



The eigenvalues of m^.ml, are to first order in a 



(,^ + 2 + ^).^ 



2 2 

p m . 



(IV.65) 



It is now easy to see the following properties of the light 
neutrinos: 

i) The neutrino mass spectrum has the inverted form. 

ii) In the no mixing limit Xg the 3 neutrino masses 
are degenerate. 

iii) The ratio between the solar and atmospheric mass 
squared differences is ~ 0.035 independent of the pa- 
rameters in good agreement with experiments. 

iv) The experimentally observed atmospheric mass 
squared difference can be used to fix the mass param- 
eter m for the light neutrinos and to estimate the mass 
parameter F'f2>3} for the heavy neutrinos: 



^{2,3} 



1 

V2x, 



0.035 eV 



1.8 • 10'^ GeV . (IV.66) 
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Here we used for Xg the estimate piL49p . Thus, as 
expected from (|II.24|1 . the largest heavy neutrino mass 
^{2,3} YiQjg indeed a value close to the intermediate mass 
scale Mi ! 

v) In the approximation used so far, and without tak- 
ing renormalization effects into account, the neutrino 
mixing matrix is of the bimaximal form. 

In the next subsection we take renormalization into 
account by running the mass matrix down from Mi to 
AIz- Moreover, the effect of xa and the correction terms 



(emerging from U) will be included. 

Neutrino masses and mixings with RG effects 

Now we make full use of Eq. pV.6ip and take the 
renormalization effects discussed in Appendix IB] into ac- 
count. Neglecting higher powers of the small quantity cr, 
the light neutrino 3x3 mass matrix at the scale Mz has 
the form 



p{l + n) 

— i 
i 



—1 1 

p{l + r2) -2xA -i^r23 + XA 

-if ?'23 + XA P 



u 



(IV.67) 



where ri,r2,r23 denote renormalization factors derived 
in Appendix |B] [see Eqs. (|B.3p . (jB.4|l ]. Their numerical 
values are 

n = 0.1325 , ra = 0.0582 , rza = 0.9426 . (IV.68) 

Without renormalization the value of p has no influ- 
ence on the neutrino mixing pattern. Now, however, p 
plays a role. It should not be large in order not to change 
the mixing pattern drastically. In fact. F^^'^^ / F^^'^^ 
must be as small as Xg in accord with p w —zjxg (|II.24p 
and our value for Xg , e.i. \p\ must be around 1. As 
a consequence, the inverted spectrum found above still 
prevails. 

With the following choice of parameters we obtain a 
very good description of the known neutrino properties. 
In this fit we adjust xa and p and decrease m slightly. 
For Xg we take the value used already. 

m = 0.0324 eV, p== 1.132, 

Xg = 0.033 , XA = -0.0601 . (IV.69) 

For the neutrino masses one gets 

mi = 0.0616 eV , ma = 0.0623 eV , 

ma = 0.0374 eV . (IV.70) 

Accordingly, the mass square differences are 

Am^^i = m^ - m^ = 8 • 10"^ eV^ , 
lAm^J =Tfil-mi = 2.5 • lO"'^ eV^ . (IV.71) 

For the mixing angles emerging from the neutrino mass 
matrix one finds - 35.3°, 6*^3 ~ 39.4° and ~ 
3.2°. By including the effect of the diagonalization of the 
charged Icpton mass matrix as obtained in section IIIIl 
the neutrino mixing angles and the CP violating phase 61 



(appearing in the neutrino oscillation amplitude) become 

012 ~ 34° , 023 ^ 43° , 013 ~ 6.3° , 

61 ~ 67° . (IV.72) 

The angles and the phase are given according to the stan- 
dard parametrization. These mixing angles together with 
the mass squared differences pV.71|l are in perfect agree- 
ment with global fits to neutrino oscillation data [l3|. 
These results are very satisfactory. For the neutrino less 
double /3-decay parameter one finds 

|(m;3^)| ~ 0.046 eV . (IV.73) 

We note that there is not much freedom within this Eg 
model to get larger or smaller values for the neutrino 
masses or the Oj^2/3-decay parameter, or to change the 
neutrino hierarchy. We already observed that F^^''^^, the 
mass of the heavy third generation right handed neutrino 
turned out to be amazingly close to Mi , the point where 
the electroweak gauge couplings of the standard model 
meet. 

As a consequence, we know now, at least approxi- 
mately, the masses and mixings of all heavy fermions. A 
direct diagonalization of the 15 x 15 neutral lepton mass 
matrix reveals in particular, that the two lightest of the 
'right handed' neutrinos (i.e. the first generation) turn 
out to have masses of only w 700 GeV. This is caused 
by the super strong hierarchy valid for these particles. 
On the other hand, the masses of the first generation of 
heavy quarks and heavy SU(2)l leptons are much larger 
(w 10® GeV). We also find, that the neutrino mixing ma- 
trix for the light neutrinos is not strictly unitary. The 
light neutrinos mix to about 2% with the neutral leptons 
L\ and Lj. The two lightest of the heavy neutrinos (a 
combination of the first generations oi i) = L2 and L3) 
mix with an amplitude of about 8 • 10"'' with the first 
generation of the light neutrinos. These heavy neutrinos 
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could be detected by their decay to the Higgs field {H")\ 
[from the SU{2)l doublet ((i?")}, (i?")2)] and a light 
neutrino provided the mass of the Higgs field {H^)l lies 
below the heavy neutrino masses. If not, they can decay 
into W and Z gauge bosons and a light lepton via the 
mixings given above. 

Our phenomenological treatment gives a clear picture 
of the spectrum and mixings of the standard model 
fermions and of their heavy Eq partners. A quantita- 
tive fit reproducing all known properties of the standard 
model charged and neutral fermions could be performed. 
Several not yet measured quantities are predicted. How- 
ever, the model is still incomplete, because an under- 
standing of the scalar sector, in particular of the vac- 
uum expectation values of the Eq Higgs fields, is not yet 
achieved. 



V. SUMMARY 

In this work we have addressed the problem of fermion 
masses and mixings. We took the gauge group Eq for 
grand unification augmented with the generation sym- 
metry SO{3)g X Pg. The fermion fields of the standard 
model are taken to be in the (27,3, representation of 
Eq X 5*0(3)3 ^ Pg- The scalar fields transforming un- 
der the generation symmetry obtain vacuum expectation 
values by a complete spontaneous symmetry breaking. 
The corresponding values provide for the hierarchy of 
the fermions, for their mixings and CP-violation. The 
generation mixings occur in conjunction with the SU (3) 
U-sj>m mixing of the standard model fermions with their 
heavy partners. To have a renormalizable model we had 
to introduce heavy Dirac fields. Integrating them out 
led to a very strong hierarchy of the right handed heavy 
neutrinos with important consequences for the light neu- 
trino properties. The onset of the intermediate symme- 
try SU{3)l X SU{3)r X SU{3)c occurs at the well-known 
meeting point of the electro-weak gauge couplings of the 
standard model. This scale determines the masses of the 
heavy neutrinos and turns out to be in full accord with 
the measured mass splittings of the light neutrinos. Also 
the formation of the Yukawa coupling matrices deter- 
mining the fermion hierarchies can happen at this scale. 
Our model needs only few fit parameters to reproduce all 
known masses and mixing properties of the fermions. Be- 
cause of the unique use of the up quark hierarchy G and 
the antisymmetric mixing matrix A for quarks, charged 
leptons and neutrinos, the mixing in the quark sector 
determines the masses, the splittings and the main part 
of the mixings in the neutrino sector. We obtain for 
the light neutrinos the inverted hierarchy and - for the 
unrenormalized case - bimaximal mixing. Taking renor- 
malization group effects and the mixings coming from 
the charged lepton sector into account and by adjusting 
the small parameter xa, the slight change to the exper- 
imentally observed mixing angles can be achieved. The 
entire spectrum of light and heavy fermions can be esti- 



mated. The two light neutrinos with 'solar mass splitting' 
have masses ~ 0.06 eV, while the lightest one weights 
~ 0.04 eV. The angle 6*13 is predicted to be « 6°. The 
CP violating phase di appearing in neutrino oscillations 
is large, around 70°. For the parameter relevant for neu- 
trino less double /3-decays we find ~ 0.046 eV. The two 
lightest members of the heavy neutrinos have a surpris- 
ingly low mass of about 700 GeV. The neutrino mixing 
matrix is not fully unitary, mixings to some high mass 
states have a magnitude of about 2 percent. 
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APPENDIX A: GAUGE COUPLING 
UNIFICATION AND RENORMALIZATION OF 
FLAVOR MATRICES 

1. Gauge Coupling Unification 

In the previous section we have fixed the field content 
of the fermion sector which are non singlets under Eq and 
contribute to the gauge coupling running. Also, scalar 
representations play an important role. In this section 
we will give the full list of fields appearing at appropri- 
ate energy scales. This allows to study gauge coupling 
unification and to examine whether or not perturbativ- 
ity is kept up to the Afpi ~ 2.4 • 10^® GeV (the reduced 
Planck mass). 

We assumed that the scales e| and £2 are both equal 
to Ml/ g2{Mi) (at Mi the couphngs gi and 32 meet). 
Below this scale the field content consists of the three 
fermion generations of the standard model together with 
two Higgs doublets. One light (up type) Higgs doublet 
hu = ((ff")}, {H^'f^) comes from 77"(3,3, 1) and has the 
VEV {{H'^)\) = e\. The second one, the down type Higgs 
doublet hd = {{H%,[H'^)l), emerges from iJ'*(3,3, 1) 
with the VEV {{H'^)l) = e\. Therefore, below the scale 
M/, the b- factors relevant for the gauge coupling con- 
stants gi, 52, 53 affected by these states are: 




(A.l) 



In addition, there are extra D and L fermion states from 
the 27" fermion representation. Their masses have the 
scales Af/((T'*, (T^, 1). The corresponding b-factors valid 
for each generation of D's and L's are 

{h^M. hf^ (^^,0, ^) , &3) - (1 1 0) . 

(A.2) 

contributing at different scales. Above Mi, all compo- 
nents of 27" can be considered light. Together with these, 
there are 'light' H"(3,3, 1) and H^(3,3, 1) scalars. The 
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b-factors corresponding to all these states are 

ibL,bR,bc)''' ^{-4,-4,-5) . (A.3) 

Moreover, at the scale M [the masses of the states 
Fa (27), (27) ] which we identify here with M/, the RG 
6-factors receive the additions bf = (12, 12, 12). Between 
the scales Mj and Mqut we have L ^ R (I?-symmetry) 
which insures the equality 5l(m) = gnifJ-)- It starts at 
Mj, the meeting point of 171 and (72- At some scale Mg 
(A// < il/g < A/gut) the scalar states (6,3,1), (3,6,1) 
[arising from iJyi(351)] have to come in. We take them 
with the common mass A/g- The leptonic states from 
Pa,F°' [Lp' + Lp,)a are also assumed to have masses 
below the GUT scale. For convenience (see the discussion 
below) we will take their masses equal to Afg- AH states 
with mass Afg are important for the successful gauge cou- 
pling unification. The corresponding b-factors are 



(&L,6H,M'=(y,y,0 



(A.4) 



By taking A/g = 3.47 • 10^^ GeV unification is completed 
at AfcuT = 2.8 • 10^"^ GeV. The mass scale A^g has been 
chosen in such a way as to insure the perturbativity of 
the unified gauge coupling acuT up to the Planck scale 
A/pi. The plot which shows the unification of the gauge 
couplings is given in Fig.[TJ It has a 'Concorde' shape 
pretty similar to the one first obtained in Ref. For 
the case considered here we have 



the scale (X33) — M after integrating out the states 
F, F. In order to have RG stability in the neutrino sec- 
tor we need to have a rather small value for the r coupling 
constant At- This is satisfied because of the small value 
of z. In a notation appropriate for the G333 symmetry, 
the Yukawa interaction involving the G matrices reads 

G^'5QLQfli?"(3, 3, 1) + zG««0LQKi^''(3, 3, 1) + 

izG^^LL7J'^(3,3,l) . (A.7) 

Since the up type quark masses are generated by e\ C i?" 
and the down type and charged lepton masses through 
e| C H'^, we have 



At = (GW)33, A, = z(GW)33, a,: 



:z(G")33 • (A.. 



We estimated zC"^ ~ 0.02 in pil.48p using ~ 76 GeV. 
Thus, RG effects for A,- and Af, are small and the only 
strong scale dependence below Mi occurs for At . For 
a detailed description we will now use the appropriate 
notation for all Yukawa couplings. Namely, e\G 
e}G™', zelG zel{G'^^,G'''^^), f^A f^A'^^ etc, 
and (G™)33 = At, z(G'''^")33 = A^, 2(G'="'=+)33 = A,, 
(A^'i)23 = iAf /2 etc. 

In the charged lepton sector we had already fixed some 
Yukawa couplings and products of Yukawa couplings with 
the corresponding VEVs. The values for At, Ah and At 
are 



Mi ~ 1.27 • 10^3 GeV , A/g ~ 3.47 • 10^^ GeV , 
A/gut ^ 2.8 • 10^^ GeV , a'lj^iMcvT) ^ 24.1 . (A.5) 

At and above A/gut we include, together with 
the states mentioned above, the remaining states of 
Fq,'(27),f" (27) in order to complete the Fg states. To- 
gether with all Fg gauge bosons, the complete Fg scalar 
multiplets H{27), H{27) and iJ^(351), also the scalar 
field i/(650) is taken into account. The role of the 650- 
multiplet is to break Fg down to G333 at the scale A/gut ■ 
The b-factor above A/qut, corresponding to all these 
states, is bpe — 63. With this input we can evolve the 
Fg unified gauge coupling above AfGUT and compute its 
value at the Planck mass: 



^GUT 



(A/pO = a^^T(A/GUT) - § In . (A.6) 



With the values of A'/gut and aGUT(^^GUT) given in 
()A.5P we find from (|A.6P aGUT(A'/pi)/47r ~ 0.03 (which 
may be viewed as an effective loop expansion parameter) . 
It is seen, that within our approach the gauge couplings 
remain perturbativc up to the Planck scale. 



2. Flavor Coupling Renormalization 



At(A/2) = 1.30, Ah(Af2) ~ 4-10"^ Xr{Mz) 2-10^^ 

.(A.9) 

The Yukawa coupling in front of the matrices A is de- 
fined according to 



A 




(A.IO) 



In our estimates we take for all Xa's the same value, 
however. 

In our model in which the mixing of Fg flavors is a U- 
spin mixing effect, one can expect that the VEVs of e|, 
/I and which appear at the weak scale have roughly 
the same magnitude. This implies using (jIII.47p and the 
discussion at the end of section [Till — gl — 76 GeV. 
Together with the fit values quoted in pll.37p . (1111.44^ 
all VEVs encountered are known at least approximately. 
Using A^^ « A^^ sa 1 and A^ = 1 we exhibit them here 
all together: 

el = 131.5 GeV, /| ~ -0.23 GeV, gi ~ -0.19 GeV, 



/I ~ 39 GeV 



9l 



76 GeV 



/I 92 



^ = -4.16- 10-6 



(A.11) 



Let us start with the renormalization of the top, bot- 
tom and tau Yukawa couplings. They are generated at 



It is important that the couplings Af, and At are small in 
this scenario. It implies that their RG dependence can be 
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ignored. We could select A^i ~ 1 because in the RG treat- 
ment at scales below Mj this coupling does not explicitly 
enter. The reason is, that the scalar states (6, 3, 1) and 
(6, 3, 1) decouple already at the higher scale Mg and thus 
cannot appear in appropriate loops. The only renormal- 
ization effects relevant below the scale Mj are due to Xt 
and the gauge couplings 51,2,3. The renormalization of 
the flavor matrices is performed in the following para- 
graph taking this fact into account. 

The flavor matrices G and A as introduced in section 
H] did not include RG effects. We start with the study 
of the renormalization of the G matrices which are gen- 
erated at the scale Mj. At this point the spacing of 
different channels (i.e. uii, dd, e~e^ etc) is universal. 
We suggested for the splitting equal spacing between 
first and second and between second and third genera- 
tions p.2p . At the scale Mz, however, for different chan- 
nels different RG factors emerge affecting the spacings. 
The G-couplings can be divided into two groups. One 
group corresponds to couplings involving i?^7. The sec- 
ond group corresponds to couplings involving Hi^j- At 

Mz they are split as follows: G^m f G"", G^^, G^^^ 



exp 



zGnd ^ z[G'^'^,G^ (we will keep only those 

quantities which are relevant for us). Our aim is to cal- 
culate the relative factors between appropriate entries in 
these matrices. The relative RG factors between the 3rd 
generation components of G"",G'^'' and G^ ^ are de- 
scribed by the coupling constants Aj, Ab, A,-. However, 
different generations have different Yukawa interactions. 
Because we can take z « 0.02 all Yukawa couplings in- 
volving this prefactor do not participating in the change 
of the spacings. With the boundary condition that all 
G-couplings at M/ have the form (|I.2p . one gets at the 
scale Mz 



QUU 



(A/z) = Diag rj-V* , r^V^ , 1 



Qdd ^ 

— (Mz) = Diag(rr'/V 



G" 



G 



DD 



Xd 
G^^ 



Xl 



-(A/z) = Diag(a^ , l) , 

(Mz) = Diag (Kfa^ , n^a^ , l) , 
(Mz) = Diag (k^ ^ ^l^2 ^ ^) _ (^^^3) 



167r2 



K2 

with t = 

Ml = 



= exp 

= exp 
In/i , 



167r2 



1 



167r2 

= 



f 2 
, (10^1 

InMz , 



9 2w, 



U = In , 
^3 = Mi. (A. 13) 



In G^^ and in G^^ the renormalization factors are due 
to gauge interactions. The reason is that the heavy 
states Moi and M^. decouple at different scales. Be- 
low the decoupling scale the mass of the corresponding 
state 'freezes' out and does not run, and the leg in the di- 
agram involving the decoupled state will not be 'dressed' 
by the loop. 

Now we turn to the RG of the A couplings: 
At low energies instead of one A matrix we are 
dealing in general with six A matrices: A — > 

\dd j^dD J^Dd j^'^ e^ A& E'^ 4-E^e+ 



The appropriate RG factors are 



scripts indicate the channels they couple to. Gompared 
to the G-couplings, the A matrices are formed at rela- 
tively high scales. These matrices are the only source for 
generation mixing, but through loop diagrams involving 
the scalar states i7^(6,3, 1) and i?^(6,3, 1), these ma- 
trices could nevertheless induce off diagonal entries into 
the G matrices. However, if the latter states are already 
decoupled when the A matrices are formed, the corre- 
sponding loops will not contribute to renormalization. 
This is one of the reasons why we take the masses of the 
Lpi + Lpi states equal to Mq. It insures that in the lep- 
ton sector non diagonal A-couplings do not influence the 
G matrices through loop effects. The A-matrices corre- 
sponding to the quark channels are formed at the GUT 
scale keeping the model intact since quarks do not couple 
with H\,H\ . 

For the renormalization of the A-matrices the relevant 
energy range is Mz — Mi where the G-matrices are al- 
ready formed (above the scale Mi the matrices A do not 
change their structural form) . Only the coupling Xt is rel- 
evant. Since only the quarks = (u, d)^ are connected 
with At, solely the (3,2) elements of the matrices A'^'^ and 
A'^^ will pick up relative RG factors. The spacing of the 
other entries and of the remaining A matrices will remain 
unchanged. Therefore, we have at the scale AIz 
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exp 



exp 
1 

16^ 
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with 
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\e-E+ 
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(T —a 
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1/2 
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-1/2 
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■^A 



(A.14) 



The numerical values of the RG factors appearing above 
are 



K 



n = 0.717 , 
° - 1.299 , 
of = 1.137 , 



KA = 1-117 , 
- 1.127 , 
Kk = 1.066 



(A.15) 



phenomenology. These fields include the G333 singlet 
5+(l,l,l) state. The VEV {S+) = Mgut/^gut gives 
the breaking Eg — > G333 in such a way that the L — R 
symmetry V is unbroken at high energies. This way a self 
consistent gauge coupling unification is guaranteed. The 
additional SM singlet VEVs are much smaller than 63 . 
They arc important for the fcrmion sector. The leading 
role for the electro-weak symmetry breaking is played by 
the VEVs el, e\, g1 and /|. The remaining electroweak 
VEVs (/I, g\) are small but still relevant for the light 
fermion masses and their mixings. 
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These coefficients are used in the analysis of the charged 
fermion mass matrices in sect. IIIII 

Finally, using the above results, we are able to fix the 
values of the SM singlet VEVs 



M 



2.27- 10" GeV 



.92 (M/) 
/I = gl = 9.5 • 10^ GcV 



(A.16) 



On the other hand, from the neutrino sector (see sect. 
IIVI) we have 



^{2,3} 



1.8 • 10" GeV 



f13'3> = x„flFl2'3} ~ 6.6 • 10" GeV 



.5 



= « -3.3-105 GeV . (A.17) 

We recall that by integrating out the TV, N states, we 
found for the mass scales in (jA.17p 



^{2,3} 



^^el4 , ^ _%^^(e3)2^ (A.I8) 

Mm Mjv ^ ' ^ ' 



(the couplings X^, Ajv come from the vertices Xj^FW N, 
XnFH^ N). The comparison of fOe)) with (|XT7l) allows 
the estimate 



el ~ 8.5-10" GeV. 



AatAjv p I 1 

(A.19) 

One recognizes again the importance of the scale il// = 
g2{Mi) 63 for the heavy fermions and in particular for 
the 'right handed' neutrinos. 

The VEV configuration discussed here is consistent 
with the charged fermion sector as well as the neutral 
lepton sector. The following picture of the symmetry 
breaking of our grand unified model emerges: In a first 
step Eq is broken by an H{650) scalar field. It reduces 
Ee to the group G333 = SU{3)l x SU{i)R x SU{3)c. 
At Mj ~ g2{Mj)el the G333 symmetry breaks down to 
the SM. The role of i?(650) fields is important for the 




We need to define rrii, at the scale Mz ■ To calcu- 
late the corresponding renormalization effects it is conve- 
nient to consider the combination {lA'^)^^mJJ~^ taking 
pV.6ip . PV.62|1 and (|IV.63|1 into account. Without RG 
effects this matrix has the form 



(B.l) 

In order to compute the RG factors we recall the origin 
of each entry in (jB.ip . Let's start with the pm-entries 
in the diagonal part. The elements (1,1), (2,2) and (3,3) 
emerge from dimension five {d = 5) operators at different 
scales (at Fl^-^lcr^, F^'^^^^a'^ and F^^.s}, respectively) by 
integrating out the appropriate states of L2 , (this can 
be seen from Fig. 2^). Above the corresponding scale 
only the Dirac Yukawa couplings run. After the forma- 
tion of the effective dimension five operators, i.e. below 
the corresponding scale, new RG factors are relevant. 
They describe the running of couplings of d = 5 oper- 
ators which are different from the. running of the square 
of the Dirac Yukawa couplings. Therefore, relative (mis- 
matching) RG factors between different entries have to 
be applied. In order to calculate these RG factors, we 
need to go through the entire energy interval consider- 
ing the appropriate thresholds. After performing these 
procedures stepwise from M/ down to Mz the diagonal 
pm-entries at Mz can be parametrized as follows 



■ r2) 



and 



pm{l + ri) , pm{l 
where the RG factors ri.2 are given by 



pm 



(B.2) 



ri = exp 

r2 = exp 
with 



1 



167r2 
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167r2 



— 9? ) dt 

+ TT^ai ) 



In Pi. 



P22 = F^^'^^a' 



3 2 
2-92 

3 2 

2^2 ■ -^Q.; 

^ j;^{2,3}^8 ^ 
M33=F^'''> 



- 1 



(B.3) 
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(when RG effects are ignored these factors are zero and 
the diagonal pm-part of the matrix (jB.ip is a unit matrix 
in flavor space). In the parametrization of (jB.2p we did 
not write any RG factor for the third term pm because 
an overaU factor can be absorbed by the parameter p (or 
the scale m) . Only the relative RG factors given in (|B.3p 
are needed. 

Considering now the remaining entries of , it is easy 
to verify that the (1,2), (1,3) elements together with the 
2mxA and tuxa entries of (2,2) and (2,3) are formed at 
the same scales after the subsequent integration of the 
heavy neutral states (this can be seen from the second 
and third terms of pV.58|) and the corresponding di- 
agram in Fig. IHb). Therefore, no relative RG factors 
between these entries will appear. The only additional 
factor occurs for the — i'Tif term in the (2,3) and (3,2) 
elements. We can parameterize this term at the scale Mz 



by — im^r23. The RG factor can be calculated from 
the equation 



exp 



(B.4) 



Taking all RG factors into account, the matrix (|B.1|) 
at Mz takes the form 

(p{^+ri) -I I \ 

-i p(l + r2)-2a;^ -if r23 + a;^ -m . (B.5) 

i ^if ^23 + XA P ) 

Finally, one obtains the neutrino mass matrix at the scale 
p = Mz as shown by Eq. pV.67p . 
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FIG. 4: Diagrams responsible for neutrino d = 5 operators. 



